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Abstract 

oc 

We consider the noncolliding Brownian motion (BM) with N particles starting from 
the eigenvalue distribution of Gaussian unitary ensemble (GUE) of N x N Hermitian 
random matrices with variance a 2 . We prove that this process is equivalent with the 
time shift t — > t + a 2 of the noncolliding BM starting from the configuration in which 
all N particles are put at the origin. In order to demonstrate nontriviality of such 
equivalence for determinantal processes, we show that, even from its special conse- 
quence, determinantal expressions are derived for the ensemble averages of products 

■ of characteristic polynomials of random matrices in GUE. Another determinantal pro- 
| cess, noncolliding squared Bessel process with index v > — 1, is also studied in parallel 

with the noncolliding BM and corresponding results for characteristic polynomials are 

■ given for random matrices in the chiral GUE as well as in the Gaussian ensembles of 
O , class C and class D. 

Keywords Characteristic polynomials of random matrices, Noncolliding diffusion 
processes, Determinantal processes, Brownian motions and squared Bessel processes 

> : 

■ 

h : 1 Introduction 

■ 

We consider iV-particle systems of the one-dimensional standard Brownian motions (BMs), 
X(t) = (Xi(t),X 2 (t), . . . , Xjy(t)),t > 0, and of the squared Bessel processes (BESQ) with 
index v > -1, X (u) {t) = (X["\t), X% \t), . . . , X%\t)), t > 0, both conditioned never to 
collide with each other, N G N = {1,2,3,...}. The former process, which is called the 
noncolliding BM [26], solves the following set of stochastic differential equations (SDEs) 

dXj(t) = dBjit) + X(t) d -X (tY 1 ^ j ^ N > (L1) 
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with independent one-dimensional standard BMs {Bj(t)}j 
noncolliding BESQ [29] , does the following set of SDEs 



N 

j=l, 



and the latter process, the 



dXf\t) = 2^Jxf\t)dB j (t) + 2(v + l)dt 




xf\t)-xt\t) 



1<J<N, t>0, (1.2) 



where {Bj{t)}^ =l are independent one-dimensional standard BMs different from {Bj{t)}f =1 
and, if — 1 < v < 0, the reflection boundary condition is assumed at the origin. (See 
[16, 10, 34, 31, 41, 48, 25, 49, 46, 36, 44, 43, 37] for related interacting particle systems.) Let 
R be the collection of all real numbers and R + = {x G R : x > 0}, and consider subsets of the 
iV-dimensional real space R^, = {x = (rr 1? a; 2 , ■ ■ ■ , x^) G R^ : X\ < X2 < • ■ ■ < %n}, and 
= {x = (xi, . . . , xn) G R+ : xi < ■ ■ ■ < xn}- The former is called the Weyl chambers 
of types Ajv_i. We can prove that, provided X(0) G and X M (0) G W£, then the 
SDEs (1.1) and (1.2) guarantee that with probability one X(t) G and X {u \t) G 
for all t > 0. In both processes, at any positive time t > there is no multiple point at 
which coincidence of particle positions Xj(t) = Xk(t) or Xj V \t) = xj?\t) for j ^ k occurs. 
It is the reason why these processes are called noncolliding diffusion processes [30]. We can 
consider them, however, starting from initial configurations with multiple points. An extreme 
example is the initial configuration in which all N particles are put at the origin. In order 
to describe configurations with multiple points we represent each particle configuration by a 
sum of delta measures in the form £(•) = Ej=i ^ 3 (')> where with given y G R, S y (-) denotes 
the delta measure such that 5 y (x) = 1 for x = y and S y (x) = for x ^ y. Note that, by this 
definition, for A C R, £(A) = f A £(dx) = Yl,i<j<N-xjeA 1 = the number of particles included 
in A. (The above mentioned example is then expressed by £(•) = N5 (-), which means that 
the origin is the multiple point with all N particles.) For a given total number of particles 
iV G N, we write the configuration spaces as VJUn = {£(•) = J2f=i ^j(') : x j e R, 1 < j < N} 
and M% = {£(•) = E^iM') : X J G R+,1 < j < N}. We consider the noncolliding BM 
and the noncolliding BESQ as 9Jtjv-valued and 971^-valued processes and write them as 



respectively. The probability law of E(t, •) starting from a fixed configuration £ G 9Jtjv is de- 
noted by P^ and that of E,^(t, ■) from £ G 931^ by P^, and the noncolliding diffusion processes 
specified by initial configurations are expressed by (E(t),t G [0, oo),P^) and (E^(t),t G 
[0, oo), P£), v > -1, respectively. We set Wl Nfi = {£ G Wl N : £({x}) < 1 for any x G R}, and 
Wljj = {£ G Wljj : £({x}) < 1 for any x G R+}, which denote collections of configurations 
without any multiple points. 

In order to dynamically simulate the random matrix ensemble called the Gaussian unitary 
ensemble (GUE), Dyson considered the NxN Hermitian matrix- valued BM and showed that 
its eigenvalue process satisfies the SDEs given by (1.1) [14]. This eigenvalue process is called 



N N 
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Dyson's BM model with parameter (3 = 2 or simply Dyson's model [47, 28]. The equivalence 
between Dyson's model and the noncolliding BM, (E(t),t G [0, oo),P ? ), implies that the 
random matrix theory [39, 15] is useful to classify and analyze noncolliding diffusion processes 
[24, 30]. In particular, if the initial configuration is given by £ = NSq, this equivalence 
concludes that, for any t > 0, the particle distribution of S(t) is equal to the eigenvalue 
distribution oi N x N random matrices in the GUE with variance t. (Note that in the usual 
GUE the mean is set to be zero.) Here the probability density function (pdf) of the GUE 
eigenvalues with variance a 2 is given by 

a ~N 2 / \ x \2\ 

/^(0 = ^^exp M^) 2 , (1-3) 



C N L V 2ff2 , 

i = X^Li fix, £ VRn, xi < x 2 < • ■ ■ < x N , where C N = (2n) N/2 Ylf =1 r (j) with the gamma 
function T(z) = J °° e~ u u z ~ 1 du, \x\ 2 = J2j =1 x 2 , and 

h N (x) = [J (xk-Xj). (1.4) 

l<j<k<N 

The expectation of a measurable function F of a random variable S G 971^ with respect to 
(1.3) is given by 

E Nta *[F(Z)] = [ F{i)n N ^{t)dx = ±- I F(^ N ,A0dx (1.5) 
Jw* Jr n 

with setting £ = J2j=i $x j ,x= {xi, ■ ■ ■ , x N ), where dx = Ylf =1 dxj. 

Assume v G N = N U {0} and let M(t) be the (N + u) x N complex matrix-valued 
BM. Then the N x N matrix-valued process L(t) = M(tyM(t),t> 0, where M(i)t denotes 
the Hermitian conjugate of M(t), is called the Laguerre process or the complex Wishart 
process [11]. The matrix L(t) is Hermitian and positive definite, and Konig and O'Connell 
proved that the eigenvalue process of L(t) satisfies the SDEs given by (1.2) [34]. Again by 
the random matrix theory [39, 15], this equivalence concludes that the particle distribution 
of (E^(t),t G [0, oo), ¥^ 5 °) at any given time t > is equal to the distribution of squares 
of singular values of (N + v) x N complex random matrices in the chiral Gaussian unitary 
ensemble (chGUE) with index v G No and variance t. Here the pdf of squares of singular 
values in chGUE with index v G N and variance a 2 is given by [51, 50] 

-2N(N+v) N 

= Hv) \\{x^^ 2 )h N {x) 2 , (i.6) 

C N j=l 

£ = Ef=i<^ eli,0<Xi<-<%, where = 2 N ^l[f=i^UW + The 
expectation of a measurable function F of a random variable S G OJt^ with respect to (1.6) 
is given by 

E&^S)] = / + F^AOdx = ± I F(0tiU0dx (1.7) 
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with setting £ = J2j=i S Xj , x = (xi, . . . , xjv). 

Let n{2N) be the space of 2N x 2iV Hermitian matrices and sp(2N,C) and so(2iV,C) 
be the spaces of 2N x 2N complex matrices representing the symplectic Lie algebra and the 
orthogonal Lie algebra, respectively. Altland and Zirnbauer introduced the Gaussian random 
matrix ensembles for the elements in U C {2N) =spn U{2N) and in U D {2N) = so nH(2N), 
which are called the Gaussian ensembles of class C and class D, respectively. The eigenvalues 
of matrices both in the class C and class D ensembles are given by iV pairs of positive and 
negative ones with the same absolute value. The pdfs of the squares of N positive eigenvalues 
are given by (1.6) with v — 1/2 for the class C ensemble and with v = —1/2 for the class 
D ensemble, when the variances are a 2 [2, 3]. If we consider the % c (2iV)-valued BM and 
% D (2iV)-valued BM, the squares of each N positive eigenvalues satisfy the SDEs (1.2) with 
v = 1/2 and v = -1/2, respectively [24]. See also [10, 31, 41, 6, 35]. 

In the present paper, we consider the noncolliding BM whose initial configuration is dis- 
tributed according to the pdf (1.3), denoted by (S(t), t G [0, oo), P Mjv - ct2 ), and the noncolliding 
BESQ starting from the distribution (1.6) with not only i/GN but with v > — 1 generally, 

denoted by (S M (t),t G [0, oo), P^V- 2 ). We prove that for any iV G N, a 2 > the following 
equalities are established, 

(S(t),tG[0,oo),P^ 2 ) ^ (H(t + a 2 ),tG[0,oo),P^>), 
(SM(t),tG[0,oo),P^ 2 ) ^ (SM(t + a 2 ),tG[0,oo),P^'), *>-l, (1-8) 

f d 

where = denotes the equivalence in finite dimensional distributions (see Theorem 2.1 and 
Remark 1). We would like to emphasize that these equalities are highly nontrivial and in or- 
der to demonstrate it we show in this paper that, even from very special consequence of (1.8), 
determinantal expressions of ensemble averages of products of characteristic polynomials of 
random matrices are derived. See [8, 17, 40, 9, 18, 19, 1, 7] for extensive study of charac- 
teristic polynomials of random matrices, especially in the connection with the Riemann zeta 
function [32, 33, 21]. 

We write the expectations of measurable function F of N x N random matrices {H} in 
the GUE, of (N + v) x iV random matrices {M} in the chGUE with v G N , of 2iV x 2iV 
random matrices {H c } in the class C, and of 2N x 2N random matrices {H D } in the class 
D as {F(H)} GVE{N ^2 h (F(M)) chGUE ( Ar ^ (T 2), (F(// c )) classC(2A r iCT 2 ) , and (-P(// D )) classD(2iV)(T 2 ) , re- 
spectively, where a 2 denote the variances of these four kinds of Gaussian ensembles. Then 
for m G N, ct G C m the ensemble averages of m-product of characteristic polynomials of 
random matrices are defined as 

I m \ I m N \ 

M GVE {m,a;N,a 2 ) = / f[ det (a n I N - H)\ = / JJ JJ(o; n - Xj) \ 

\n=l I GUE(iV ,<r 2 ) \n=li=l / GUE (jv j<7 2) 

(1.9) 
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M^ VE (m,a;N,a 2 ) = ( J] det(a n I N - M*M)^ 



\n=l 
m N 

nn< 

\n=l j=l 



a„ — Ka 



chGUE(Af,y,cr2) 

, V G N , 



(1.10) 



chGUE(7V,^<r 2 ) 



and for Jj = C and D 



M classfl (m, a, 2N, a 2 ) = ( JJ det(a n / 2JV - # tt ) \ 

/ class 



Vn=l 



m N 

\n=lj=l I c lasst(2N,a 2 ) 



classft(2iV,(T 2 ) 
m N \ 

nii(^-s 2 )) .a-") 

n=lj=l I classtt(2Ar,<T 2 ) 



where U denotes the t x t unit matrix, (Ai, . . . , Ajv) are the eigenvalues of H, {k\ , . . . , /c^-) are 
the eigenvalues of M^M, (ei, . . . ,en, — £i, . . . , — £at) are the eigenvalues forming "particle- 
hole pairing" of H^, (j = C or D. For (1.10) remark that each M in chGUE has such a singular 
value decomposition that M = WKV, where U G U(iV + u), V G U(iV), 



AT 



with K = diagjfti, k 2 , • ■ ■ , k n }, (k±, . . . , k n ) G W 



N 



and O is the u x N zero matrix. The diagonal elements (k±, . . . , kjv) of K are called the 
singular values of rectangular matrix M. Since M^M = V^K T KV, the eigenvalues of M^M 
are squares of singular values (k\, . . . , k 2 n ) [20]. 
We use the convention such that 



[T/(^) = exp<^ / i(dx)\ogf{x) 



n w ({x}) 

xGsupp £ 



for £ G OJIat and a function / on R, where supp £ = G R : £({:r}) > 0}. Then (1.9) is 
given by 



M GUE (m, a; N,a 2 ) = 
with (1.5). And if we define 



niiK-x) 



.n=i xes 



;i.i2) 



n n k - x ) 



M^(m,a;iV,a 2 )=E( Vi(T2 

_n=l xes 

with (1.7) for z/ > —1, (1.10) and (1.11) are given as 

^gueK «; a 2 ) = MM(m, <*; TV, a 2 ), z/ G N , 
M classC (m, a; 27V, a 2 ) = M^ 2 \m, a< 2 >; iV, a 2 ), 
M classD (m, ct, 2N, a 2 ) = M^ 2 \m, a< 2 >; TV, a 2 ), 



(1.13) 



1.14) 
(1.15) 

;i.i6) 
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where for a = (a u ...,a m )e C m , ct^ = (aj, . . . , a 2 m ) G C m . 

We will show that, from the equalities (1.8), two sets of determinantal expressions are 
derived for (1.12), (1.14)-(1.16). 



2 Preliminaries and Main Results 

We define 

p(t,y\x) = 



exp 



2t 



t Gl\{0}, 
t = 0, 



(2.1) 



y/2ir\t\ 
5(y-x), 

for x, y G C [27]. For z G C, v > — 1, we define z v to be exp {y log z), where the argument of 
z is given its principal value; 

z v = exp |Vj log \z\ + v^Targ(2;)|j , — n < arg(^) < n. 

For v > —1, y G C, we set 



V {u \t,y\x) = I 



2 t \x 



yy/2 



exp 



(2^1)^(^ + 1) 
5(y-x), 



x + y 
2t 

,-y/2t 



t Gl\{0},ieC\{0}, 

* 6l\{0},i = 0, 
t = 0,x G C, 



where I v (z) is the modified Bessel function of the first kind defined by [52, 4] 



(2.2) 



1 z ' 

h{z) = ^ T(n + l)T(n + l + u) (2) 

n=0 



The following equality holds, 

dyp(s, z\y)p(t, y\x) = p(s + t, z\x), 
dyp( u \s, z\y)p( u \t, y\x) = p ( - u \s + t, z\x) 



J 

Jr 

/ 



(2.3) 
(2.4) 



for s,t > 0,x,z G C. The former is confirmed just performing the Gaussian integral and 
the latter is proved by using Weber's second exponential integral of the Bessel functions 
[52] with appropriate analytic continuation [29]. When t > 0,x, y G R, p(t, y\x) gives the 
transition probability density of the one-dimensional standard BM from x to y during time 
period t, and when t > 0,x, y G K+, p^\t,y\x) gives that of the BESQ with index v > — 1 
from x to y during time period i (if — 1 < 1/ < 0, a reflection wall is put at the origin). 
The equalities (2.3) with i,z6R and (2.4) with x, z G IR+ are the Chapman-Kolmogorov 
equations of these one-dimensional diffusion processes [45]. 
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We introduce the Karlin-McGregor determinants 

f(t,y\x) = det A Mt,yj\x k )], x,yeW^, 
f^(t,y\x) = det y»\t, yj \x k )}, x,yeW+, v > -1, (2.5) 

t > 0. 

For given f G Wl N , N G N, we write £(•) = £f =1 $y (•) with en < a 2 < • • • < a N . Then 
we set 

n 

6(0 = and £„(0=J> a ,(0, l<n<iV. 

j'=i 

For i > 0, x G C, we define [5, 28] 

<j>P(t,x;S) = ^-l dsp(t,x\s)- -, (2.6) 

^(^x-Z) = [ dsp(-t,is\x) T[(is-a), £eWl N , (2.7) 



aei. 



and [12, 29] 



^' + \t,x;0 = ^<f d8pM(t,x\s)= r, (2.8) 

<^(M;0 = / dspM(-Ml*)II( a - a )' ( 2 - 9 ) 

n = 0, 1, . . . , N — 1, where i = v^T, for £ G Tie, ( G N, C(C) denotes a closed contour on 
the complex plane C encircling the points in supp ( on the real line R once in the positive 
direction, and R_ = {x G R : x < 0}. And for t > 0, let 

MJ (*,!/;£)= del ^ [^3 (*,!/*; 01> ee^, 

±<j,k<N J 

^' ±) (t,y;0 = 1 del ; [^f(t,y fc ;0], !/ e W+, £ e SW+. (2.10) 

l<3,k<N J 

Since </>i ^(t, x;£) and 0^' ^,x;£) are monic polynomials of x of degree n, 

h { -\t,y;Z) = h%-\t,y;Z) = h N (y), (2.11) 



which are independent of t and £. On the other hand, the following equalities are proved 

Jj = l Xj 



(Lemma 3.1 in [28] and Lemma 3.4 in [29]). For any t > 0, y G W^,£ = EL ^ G % 



with xi < x 2 < • • • < xat, 

fMx) -h%\t,y;0, (2.12) 



h N (x) 
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and for any t > 0, y G W£, £ = Ef=i ^ e 9tt£ with < 



/ (1 %?/l*) 



? , 1 =h% + \t,v,£), (2.13) 
Iin{x) 

where, if some of the x/s coincide, the LHS of (2.12) and (2.13) are interpreted by using 
l'Hopital's rule. 

For any MeN and any increasing time-sequence < t\ < ■ ■ ■ < tu < oo, the multitime 
probability density of (E(t),t G [0, oo),P ? ) is given by 

p^i,£ (1) ;---;*m,£ (m) ) 

M-1 

= h { ~\t M ,x {M) ;£) J] /(Wi - *m;« (m+1) l* (m) )fc (+) (*i,e (1) ;0 (2-14) 

rra=l 

with £ = Ejli ^ e 9^v> a i < a 2 < • ' • < a N for the initial configuration and = 

Yl7=i ^ ( m > e 9#ivo, a;( ' m ' ) = (^i m \ • • • , ^aT'*) ^ for configurations at times t m , 1 < m < M 

x j 

[28]. Similarly, the multitime probability density of (E^(t),t G [0, oo),P|) is given by 

p^i,£ (1) ;...;*m,£ (m) ) 

M-1 

= h^-\t M ,e M) ;Z) J] / M (Wi - *m;« (m+1) |aj (m) )fc (,/ ' +) (*i,e (1) ;0 (2.15) 

m=l 

with £ = J2f=i £ 9^jv; < ai < a 2 < ■ ■ ■ < a N for the initial configuration and £ (m ) = 
Y^j=i $ x ( m ) e ^jv,0' a;<m '' ) = (^i" 1 ^; •-. 5 e for configurations at times t m , 1 < m < M 
[29]. By definitions (2.6) and (2.7), we can see that 

6 {+) (t x-N5n) = r {n+1)/2 2 ^ H (^—\e~ x2/2t 

<P ( -\t,x;N5 ) = t n/2 2~ n / 2 H n (^=y < n < N - 1, (2.16) 
where H n (x) is the Hermite polynomial of degree n, 

,k (2x) 



[n/2] , 

' j kUn 

k=0 



*„(*) = ^(-l) B(| . a)| 

, r f ~y 2 /2 ^ v \ r e 2zx-z 2 

= T' 2 dy e —=(iy + V2x) n = ^ I dz 6 

JR V27T 27T1 J C{S()) 



~n+l 

C(5 ) 



(2.17) 

where [r] denotes the largest integer that is not greater than r G R + , and by definitions (2.8) 
and (2.9), 

n-(n+u+l) . „ x v / t \ 

tiTKWV = r<»-)(-i)» F ^ T - TTy (£) e-^(|), 

4"'-»(«,i;JV*) = f"(-l)"2"»!^ (|) , (2.18) 



where Ln\x) is the Laguerre polynomial of degree n with index v, 

L»(x) = Vf T(n + u+l)x k 

n[ ' ^ [ ' r(k + u + l)(n-k)\k\ 



2m J c 

Then we can prove the following. 



- 1 L^ri du . (2 . 19) 



C(5 ) 



U" 



Theorem 2.1 For any N, M G N 7 any increasing time-sequence < ii < ■ • • < t M < oo, 
and any a 2 > 0, 

\p s (h, £ {1) ; • • • ; t M , £ {M) )] = P m ° (h + o\ . . . ; t M + a 2 , ^ M) ) (2.20) 
with £ (m) G Wl Nfi , \<m<M, and 

■ ■ ■ ; t M ^ M) )} = P ? So (h + <t 2 ,£ (1) ; ...;t M + a 2 ,e (M) ) (2.21) 

with ^ G Wl+ fi , 1 < m < M. 

Remark 1. When M-time probability density of a process is given for any M G N and 
any time sequence < t\ < • • • < tM, it is said that the finite dimensional distributions 
of the process is determined [45]. Eq.(2.20) (resp. Eq.(2.21)) means that the processes 
(E(t),t G [0,oo),P^ 2 ) and (E(t + a 2 ),t G [0, oo), W> NS °) (resp. (E^(t),t G [0, oo), P^ ct2 ) 
and (S^(i + cr 2 ),t G [0, oo), F^ s °)) are equivalent in finite dimensional distributions for any 
a 2 > 0, which is the fact expressed by (1.8). 

For x^ = {x ( ^\...,xP) G or x^ G W+ with ^ m \-) = £L 5 (m) (-) and 

N' G {1,2,..., N}, we put = (x^ , . . . , x^) G or G W+„ 1 < m < M. For 
a sequence (iV m )£f =1 of positive integers less than or equal to N, we define the (N±, . . . , JVm)- 
multitime correlation functions of (S(t),t G [0, oo),P^) and (S^(t),t G [0,oo),Pf) by 

oHt, x {1) - x (M) ) 

. M N M 

- / n „ R _ n n ■ ■ ■ ; ? (M> ) n 

J lim=i M m m=l j=jV m +l m=l V my 



M N M 



n n ^%^ (i) ;---;^an7Fivv' ^ 

T-rM K JV-iV m XX XX - LJ - A' — iV m ) 

llm=l M + m=l j=JV m +l m=l v my 

respectively. In the previous papers we have shown that for any fixed initial configuration 
£ the noncolliding BM and the noncolliding BESQ with finite numbers of particles are 
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determinantal processes in the sense that any multitime correlation function is given by 
a determinant [28, 29] 

. . .;t M ,x { ^2) = det , [K ? (t m , xj m) ; t n , x^)} 

i-_^J _ 1 * m i J-S: f^"Si 

l<m,n<M 

pi(t 1 ,x^,...-t M ,x ( i2)= det [Kl(t m , Xj^-jtn, x^)], u>-l. (2.23) 

l<m,n<M 

Here the correlation kernels are given by 

AT-l 
n=0 

Z(dx') / dy'p(s,x|xO$(£;^iy>(-t,iy'|y) 

-l(s>t)p(s-t,x\y), £e 971^,0, (2.24) 

AT-l 

K^x;^) = £ <^ + )(s, x; 0$'-\t, y; ~ l(s > t)p^(s - t,x\y) 

i{dx') [ dy'p^\s,x\x')^x',y')p^\-t,y'\y) 
-\{s>t)p { ?\s-t,x\y), £e0tt+ u>-l, (2.25) 



/ 

Jr 



where 



= J] -, (2.26) 

aee,f\{x} c 

and l(ou) — 1 if ui is satisfied and l(ou) = otherwise (Proposition 2.1 in [28] and Theorem 2.1 
in [29]). The function a;, z) is an entire function of z G C expressed by the Weierstrass 
canonical product with genus 0, whose zeros are given by supp £ fl {x} c [38, 27]. 

As direct consequences of Theorem 2.1, we have the following equalities for multitime 
correlation functions; 

E N ,Ap E (h, t M , x%>)] = P ma {h + a 2 , a#>; ...;t M + a 2 , x%>) (2.27) 

with x { ™l e W^ m , 1 < m < M, and 

e^[^(*l = + ^ • • • ; *m + °\ ( 2 -28) 

with e W^ m , 1 < m < M, for any M G N, < ti < • • • < t M < oo. Therefore, infinite 
systems of equalities between determinants of correlation kernels are obtained as a corollary 
of Theorem 2.1. 
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Corollary 2.2 For any N, M e N, < h < ■ ■ ■ < t M < oo,N m e {1, 2, . . . , N}, 1 < m < 

M, a 2 > 0, 



K . <jV det [K s (t m ,xJ m) ;t n ,4 n) )] 

l<3<N m ,l<k<N n 
l<m,n<M 

Km,n<M 



(2.29) 



E 



iV,a 2 



l<m,ri<M 

i<j<A^m,i<fe<Ar„ 

l<m,n<M 



(2.30) 



In particular, as M — 1, Ni — L < N, 



E 



del ; [K s (t,Xj-;t,a; fc )] 

l<j,ft<L 

, del ; [K^(t,Xj;t,x fc )] 



= det [K w5 °(t + a 2 ,x,-;t + o- 2 ,x fc )], (2.31) 

l<j,k<L 

= det [Kf^ + ^^ + cr 2 ,^)] (2.32) 



/io/d /or any t > 0, a 2 > 0. 



The proof of Theorem 2.1 is given in Sect. 3.1. 

The main purpose of the present paper is to show that the equalities in Corollary 2.2 are 
nontrivial even in the special cases given by (2.31) and (2.32), and from them the determi- 
nantal expressions for the ensemble averages of 2n-products of characteristic polynomials of 
random matrices are derived for any nGN. We show two sets of determinantal expressions. 
The first one is given by the following theorem. 

Theorem 2.3 For any N, n e N, a = (cci, a 2 , ■ ■ ■ , a 2n ) e C 2n 7 a 2 > 0, 

T n(2N+n) 



M GUE (2n,a;iV, a 2 ) 



x det 

l<j,k<n 



1 



QLj O n _|_fc 



h n (a l , • • • , a n )h n (a n+1 , ■■■ , a 2n ) 

H N+n (a j /\/2a 1 ) H N+n (a n+k / V^a 2 ) 



H^n^aj/V^a 2 ) H^n^an+k/V^a 2 ) 



with 



9 -n(2AT+2n-l)/2 A (N + U - £)\ 

lN ' 2n= l\(N + n-iy: 



(2.33) 



(2.34) 
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and for v > — 1 

M^(2n,a;N, a 2 ) 



l { An^ 2 T {2N+n) 



h n (ai, • • • , a n )h n (a n+1 , ■■■ , a 2n ) 



x det 

l<j,k<n 



(Xj - a n+k 



J N+n K^J 



aj/2 



(T 



L v N+n (a n+k /2a 2 ) 



L u N+n -i (%/2^ 2 ) K +n ^ (a n+k /2o- 2 



with 



/at I \| \ n-1 n-1 n+1 

^» = ( rww'J II r(w + - + 1) II (Jv + m - 1)!. 

' '' l=\ m=l 



(2.35) 



(2.36) 



By setting v G N , ^ = 1/2 and z/ = —1/2 m (2.35), the determinantal expressions are given 

chGUE' 



/° rM chGUE' M ctac and M classB through (1.1 4)- (1.1 6). 



Proof is given in Sect. 3.2. The above expressions can be simplified by using the following 
identity, which was recently given by Ishikawa et al. [22]. For n > 2, x = (xi, . . . , x n ), y = 
(yi, . . . , y n ), a = (ai, . . . , a n ), 6 = (6i, ...,b n )e C n , 



det 

l<j,k<n 



Vk - Xj 



1 aj 
1 h 





i 


Xi ■ 


x 1 


ax 


a±Xi ■ 


■ ■ CL-yX™ 




i 


x 2 ■ 


T "-l 
x 2 


a 2 


a 2 x 2 ■ 


■ ■ a 2 X2 


(_l)n(n-l)/2 


i 


Xn 


~n— 1 


a n 


Q"n,Xn 




=1 mute* 


i 


Vi ■ 


■ vr l 




hyi 


■ ■ hy? 




i 


2/2 • 


■ vt 1 


b 2 


b 2 y 2 


••w 




i 


Vn ■ 


■ Vn' 1 


b n 


bnVn 


• • b n y n n 



(2.37) 



Remark 2. For n G N,p, g G N satisfying p + q = n, and x, a £ C n , denote by V p ' 9 (a;; a) 
the n x n matrix with j-th row 



(1, Xj, ■ ■ ■ ,Xj , %, &jXj, , OjXj ). 

If g = 0, then p = n and V™'°(cc; a) = V n '°(£c) = (^ _1 )i<j,fc<n is the Vandermonde matrix 
and its determinant det K™'°(a;) is equal to the product of differences of n variables, h n (x), 
given by (1.4). As a generalization of the Cauchy determinant 

h n {x)h n (y) 



det 

i<j,k<n \ Xj +y k J n"=i IlLi^i + Vk) ' 

x,y £ C n , Ishikawa e£ a/. [22] proved the following equalities involving the generalized 
Vandermonde determinants V p ' g . Let n G N, p, q G N . For x = (xi, • • • ,x n ),y = 
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(yi, • • • , y n ), a = (a u • • • , a n ), b = (b ± , ■ ■ ■ , b n ) G C n and z = (z 1} ■■■ , z p+q ), C = (ci, ■■■ , c p+q ) G 



det 

l<j,k<n 



Vk - Xj 

_-Qn(n-l)/2 



dety p+w (^,y fe ,z;a j ,fe fe ,C) 



n?=i Ul=i(yk-Xj) 



det V M (z, CY 1 - 1 det V n+p ' n+q (x, y, z; a, b, C). (2.38) 



When p = q = 0, det V°'°(z, C) = 1 and 



1 CLj 

l h 



detV ' {xj,y k ;aj,b k ) = 

Then as a special case of (2.38), (2.37) is obtained. 
For a 2 > define 

H t (a; a 2 ) = {^j H t (-j== ] , «6l, 
L?(a; a 2 ) = (-2a 2 )^ (^L) , a G 



(2.39) 
(2.40) 



£ G N , which are both monic polynomials of a with order i. By using the identity (2.37) 
and recurrence relations of Hermite polynomials and Laguerre polynomials, we can prove 
the following second set of determinantal expressions. 



Theorem 2.4 For any N, n G N, a = (a 1: a 2 , • ■ • , a 2n ) G C 2n , a 2 > 0, 

1 



M GVE (2n,cx;N, a 2 ) 
M^l VE (2n, a ;N,a 2 ) = 

M classC (2n,a;2iV,a 2 ) = 



det 

h 2 n{Ot) l<j,k<2n 
1 



(2.41) 



det 



^2n(a) l<i^<2n 



L^.^a 2 )] , z/GNq, (2.42) 



1 



det 



1 

p2n 



rl/2 / 2 2\ 



det 



#2iv+2,-i(a fe ;a 2 ) , (2.43) 



and 



M classD (2n,a;2iV, a 2 ) 



1 



det 



h2n{.OL {2) ) l<J,fc<2n 
1 



9-1/2 / 2 2* 



h2n{OL^>) l<^fc<2n 



(a fc ; ^ 2 )] • (2.44) 
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The determinantal expressions (2.41)-(2.44) can be obtained from the general formula 
given by Brezin and Hikami as Eq.(14) in [8]. (See also [17, 40, 18, 19] and Sect. 22.4 in [39].) 
Since our new expressions (2.33) and (2.35) are independently derived in this paper, if we 
combine the present result and that of Brezin and Hikami [8], the special case of identity 
(2.37) of Ishikawa et al. [22] is concluded. 

Remark 3. In this paper, we derive the ensemble averages of products of characteristic 
polynomials of random matrices from the equalities (2.31) and (2.32). As stated in Corollary 
2.2, these equalities are special cases with M — 1 of the systems of equalities (2.29) and 
(2.30). It will be an interesting problem to clarify all the information involved in (2.29) and 
(2.30) (see [13]). 



3 Proofs of Theorems 



3.1 Proof of Theorem 2.1 

By (2.16) and the fact that H n (x/2) is a monic polynomial of order n G N and thus 
deti<j t k<N[Hj_i(xk/2)] = hj^(x), we obtain the equality 

h^\a 3 ,x;N5 ) = ^^-, ^f^G%xeWi (3.1) 

Combining this with (2.12) gives the equality 

f(t u x^\x)h^\a 2 : x;NS ) = htffa, *«; £ W(0 (3.2) 
for any fj. > 0, £ = J2f=i ^ e %K N , x^ G Wft. Then (2.14) with f = N5 gives 

^(a 2 ,e;ti + ^,e (1) ;...;tM + ^,e (M) ) 

M-1 

= h { -\t M + a 2 ,x^;N5o) ]J f^-t^x^lx^f^x^lx^^cr^x^So) 

m=l 

M-1 

= h { -\t M ,xM-0 J] /(Wi-t m ;^ (m+1) k w )4 +) (ti,( (1) ;^v(0 

m=l 

= mAOpHti,C {1) ;---;t M ,C {M) ), (3.3) 

where (2.11) was used. Integrating the both sides with respect to £ over VJUn according to 
(1.5), (2.20) is obtained. Similarly, the equality 

f^( tl ,x^\x)h^ + \a 2 ,x;NS ) = h% + \t u 0/<U0 (3-4) 

is established for any t± > 0, £ = Zljli ^ e £^ e W]y-, and then we have 

^ e (1) 5 • • • ; t M , £ (M) ) = pf V, £; f i + ^, . . . ; t M + a\ (3.5) 

Integrating the both sides with respect to £ over 9Jt^ according to (1.7), (2.21) is obtained. 
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3.2 Proof of Theorem 2.3 

First we derive the expression (2.33) from (2.31). We observe that the LHS of (2.31) is 
written as follows. 

Lemma 3.1 For any a 2 > 0, t > 0, 1 < L < N, x L = (x u ...,x L )e W£, 



E 



N,a 2 



, det i \K a (t,Xj]t,Xk)] 

l<3,k<L 



-L(2N-L) 



du l dwe \ u \ 2 l 2(y2 TT p(t xAuj)p(— t, iwAxj) 

xh L (u)h L (iw)M G u E (2L, (u ± , ...,u L , iw u iw L ); N - L, a 2 ). (3.6) 



Proof Assume that £ = J2f=i ^ w ith a = (°i> • • • > a Jv) e W^-. Then 
p(t,x L ) = ^det^jK^^tjXj-jtjXk)] 



= det 



9=1 



y]p(t,gj|a g ) / rfwp(-t,iw|x fc ) ]} 



iw — at 



i<e<N,e^ q a i ae 



l<qi<-<q N <N 



3=1 



1<tj<N,tj^qj q i 



a qi - ae, i<j,k<L 



det Mt,Xj\a qk )], 



where multilinearlity and antisymmetry property of determinant have been used. Let I at = 
{1, 2, . . . , N}, I L — {1, 2, . . . , L}. For a given ordered set of indices = {qi, (fe, ■ ■ ■ , <?l}, 1 < 
gi < • • • < < iV, we see 



h N (a) 2 n n 



iWj - 



x ii n ( a *k - a 3)( iw k - * n n 



We also note that 



r -N 2 



— (N—L) 2 



X 



T -L(2N-N') 



N\C N (N-L)\C N _ L (2tt) l /2 (N-n)\ 
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Then the LHS of (3.6) becomes 

-TV 2 



\ da^—e-\ a \ 2 ^h N {afp^t )XL ) 

J R N IWlsN 



.l)L{L- 



l)/2 a -L(2N-L) 



n L Li n 

L\ J2 / dwJlpi-t^lxj)]] / da qk e-<^ 2 

^ ^„ ^ A tJK L „_i K_ n JR 



x n n ~ o K^i^^' x ;K)] E iv-L, CT 2 

(_l)L(L-l)/2 (T -L(2iV-L) - L ; 



n n K fe -ai)^-^) 

fcei L jei L \qi, 



L / d«,IJp(-*,i^)/ ^e^l 2 /- 2 JJ H (i Wk - v .) 

x det [p(t,Xj> fc )]M GUE (2L, (t^, . . . , iu^, . . . , iw L ); iV - L, a 2 ), 
where we have replaced the integral variables (a gi , . . . , a gi ) by (t>i, . . . , -ul) = v. By definition 

" res L j=l 

where Sl denotes the collection of all permutation of (1,2, ...,L). For each r e <Sl, set 
f r (j) = Uj, 1 < j < L, that is, Vj = u T -i^, 1 < j < L. Then the above equals 

(_l)L{L-l)/2 a -L{2N-L) r r Jl. 

fo ' l2T1 L , AT 77 du dwe-M / 2CT [[pfaxjluM-t^Wilxj) 

(27r) L / 2 [\n=l( N ~~ n )- J RL J ^ L j=l 
L 

X S g n ( T ) II II ( [Wk ~ U r~ 1 {3)) 

reS L k=l \<i<Lj+k 

xM GVE (2L, (u T -i(i), . . . , « T -i(L), iwi, • • • , iw L ); N - L, a 2 ). 

By definition (1.9), Mgue^, N, a 2 ) is symmetric in a, and thus 

M GU e(2L, (u T -i(i), . . . , « T -i(L), • • • , iwi,); N - L, a 2 ) 
= M GUE (2L, . . . , u L , iw 1 , . . . , \w L )\ N - L, a 2 ). 

We can confirm that for any n > 2, x, y £ M 71 

n 

8gn(«) II II ( x * " y*U)) = (-l) [n/2] h n (x)h n (y). 

K&Sn k=l l<j<n,jj^k 

We can see that [L/2] + L(L — l)/2 is even for any L £ N. Then (3.6) is obtained. | 
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Then we consider the RHS of (2.31). By the first expression of (2.24) and the fact (2.16), 
we have the extended Hermite kernel [42, 48, 26] 

K N6 ° (t + a 2 , x;t + a 2 ,y) = ^= V — H n \ H n \ - I 

(3.7) 

t > 0, a 2 > 0, x, y G R. By the extended version of Chapman-Kolmogorov equation (2.3), 
the following integral formulas are derived. 

Lemma 3.2 For t > 0, a 2 > 0, n G N , x, y G K 

/ \ / 2 \ ("+l)/2 / \ 

■M*,*)* (^) = (^) *. 

(3-8) 

H n P("*. = ( #n I /7T7T~ oT I , (3-9) 



K^>(t + a 2 ,:r;t + a 2 ,f)= 1 /" d« / dv e^ 2 ^ 2 p(t, x\u)p(-t,w\y) 

2 JV (7V - 1JV7T JR Jr 

(3.10) 



1 

x 

u — if 

Proof Since the equality 



H N {u/V2o^) H N {iv/V2o^) 



I dup(t, x\u)p(a 2 , u\s) = p(t + a 2 , x\s) 
Jr 

holds for t > 0, a 2 > 0, x G R, s G C, (2.6) gives 

/ dup(t, x\u)^\a 2 , u- = 0i +) (t + a 2 , x- 0, 
./ii 

n G N ,t > 0,o- 2 > 0,x G M. If we set £ = iV<5 and use the fact (2.16), we obtain (3.8). 
Similarly, by the equality 



dvp{— <r 2 , is | if )p{— t, if \x) = / dv p(a 2 , s\v )p(t, v | — ix) 

Jr 

= pit + a 2 , s| — ix) = p(—(t + (T 2 ), is|x), 
i>0,(T 2 >0 ) i6l 1 and by the definition (2.7), we see 



/ dv (p { -\a 2 , if; £)p(-t, iv\x) = <P { ~\t + a 2 , x- £). 
Jr 
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If we set £ = N5 and use the fact (2.16), we obtain (3.9). Inserting (3.8) and (3.9) into the 
RHS of (3.7), we have 

K NS °(t + a 2 ,x;t + a 2 ,t) = 1 I du I dv e~ u2/2a2 p(t, x\u)p(-t, iv\y) 

v2ira 2 Jr Jr 



N-l 



x 



U 



H„ 



Here we use the Christoffel-Darboux formula (see, for example, [4]), 



if 



jv-i 1 1 
^ 2^J Hn ^ Hn ^ = 2 N (N-1)\( 



H N (x) H N (y) 
H N _i(x) H N _i(y) 



ly.(x-y) 
and then (3.10) is obtained. | 

The RHS of (2.31) is thus written as follows. 

Lemma 3.3 For any a 2 > 0, t > 0, 1 < L < N, x L = (x b ...,x L )e W£, 
det [K N5o (t + G- 2 ,Xj-,t + o- 2 ,x k )} 

l<j,k<L 

— du dw e~\ u \ 2 / 2a2 Y\ p(t,Xj\uj)p(—t,iWj\xj) 
Jr l Jr l fj[ 



1 



Het 

(2 N (N - l)\y/F) L l<j,k<L 



Uj - \w k 



H N ( Uj /V2o^) H N (iw k /V2o^) 



(3.11) 



Proof By definition of determinant, the LHS of (3.11) is given by 



(2 JV (iV-l)! v ^F) L 



^ sgn(r) / du dv< 

t( z Sl JR L J^- L 



-\U\ 2 /2a 2 



1 



Note that for each permutation r e S L 

L 



H N ( Uj /V2o^) H N (iVj/VM) 



p{-t,Wj\x T{j) ) \ . 



3=1 3=1 

We set v T -i(j) = Wj, 1 < j < N, that is, Vj = w t q), 1 < j < N. Then the above equals 

2^«iv -i)o^/2 J RL du J RL dw e ~ lul2/2a2 nw. x ^M-t, i^-)} 



x s s n ( r ) n 

reS L 3 = 1 



Uj — 1W 



H N { Uj /y/2o^) H N (iw TU) /y/2o*) 
HN.^Uj/VM) H N ^{iw T{j) /V2o^) 
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which is the RHS of (3.11). i 



Proof of Eq. (2.33) in Theorem 2.3 Since the equality (2.31) in Corollary 2.2, holds for any 
t > 0,a 2 > 0, Lemmas 3.1 and 3.3, imply that the integrand of multiple Gaussian integral 
in the RHS of (3.6) is equal to that in the RHS of (3.11). By replacing L and N — L by n 
and N, respectively, and (ui, . . . , ul, iu>i, • • • , iwl) by («i, . . . , a n , a n+ i, . . . , a2 n ) = a G C 2n , 
we obtain (2.33). ■ 

Next we prove (2.35). The LHS of (2.32) is written as follows. 
Lemma 3.4 For any a 2 > 0, t > 0, 1 < L < N, x L = (x 1 , ...,x L ) G W£, 



E 



N,a 2 



, det ; \ti£(t,Xj',t,x k )] 



= 4 du L dw u^^^'^-^M^ m +v+1 . e)iN _ 01 

x/n(u)/i t (w)MM(2L, (ui, . . . , « L , iui, . . . , iu L ); N - L, a 2 ). (3.12) 



Since we can prove this lemma in the similar way to Lemma 3.1, we omit the proof. 

Then we consider the RHS of (2.32). By the first expression of (2.25) and the fact (2.18), 
we have the extended Laguerre kernel [15] 



K™ (t + a 2 ,x;t + a 2 ,y) 



x u e -x/2(t+a 2 ) N ~! 



71! 



LI 



{2(t + a 2 )}^ ^T(n + v+l) n \2(t + a 2 )J n \2(t + a 2 ) 



x 



, (3.13) 



i>0,ff 2 >0,3;,i/6 R + . By the extended version of Chapman-Kolmogorov equation (2.4), 
the following integral formulas are derived. 

Lemma 3.5 For t > 0, a 2 > 0, n G N , x, y G R+ 



u 

" v 2^ 



a 



n+u+l 



LI 



x 



n \2(t + a 2 ) 



x u e -x/2{t+a 2 ) 



(3.14) 
(3.15) 



K™°(t + o- 2 ,x;t + cr 2 ,y) = 



(2a 2 yr(N + v) J R+ 

1 



X- 



u — v 



— I du I dvu v e^ l ^p {v \t,x\u)p ( r\-t,v\y) 

(3.16) 



L v N {uj2a 2 ) L» N (v/2a 2 ) 
L^u/20- 2 ) L^v/20- 2 ) 
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Proof Since the equality 

/ dup {u) (t,x\u)p {u \a 2 ,u\s)=p M (t + a 2 ,x\s) 

holds for t > 0,a 2 > 0,x eR + ,s eC, (2.8) gives 

dup^(t, x\u)^ + \a 2 , u; £) = ^ + \t + a 2 , x; £), 



i? 4 



n G N ,t > 0,<r 2 > 0,x G R+. If we set £ = iV<5o and use (2.18), we obtain (3.14). Similarly, 
by the equality 

/ dvp^\— a 2 , — w\v)p^ v \— t, v\x) — / dup( u \a 2 , w\u)p^\t, u\ — x) 
Jk- Jr + 

= p {v \t + a 2 ,w\ -x) =p {v \-{t + a 2 ),-w\x), 
t > 0, a 2 > 0, x, w G R+, and by the definition (2.9), we see 



If we set £ = N5o and use (2.18), we obtain (3.15). Inserting (3.14) and (3.15) into the RHS 
of (3.13), we have 

' 1 du f dvu u e- u/2a2 p {u \t,x\u)p^\-t,v\y) 

</R_ 



K™°(t + a 2 ,x;t + a 2 ,t) = 



(2a 2 Y +1 

N-l 

*E 

n=0 



T(n + u + l) n \2a 2 



u 



Here we use the Christoffel-Darboux formula [4] 

N-l 



E 



ni 



Then (3.16) is obtained. | 



r(N + v)(x-y) 



L N {x) V N (y) 



The RHS of (2.32) is thus written as follows. 

Lemma 3.6 For any a 2 > 0, t > 0, 1 < L < N, x L = (x 1 , ...,x L )e 
det [K^°(t + a 2 ,x i; t + cr 2 ,:r fc )] 

= [ du [ dw\]u v j e- u i' 2 ° 2 p^\t,x j \u j )pW{-t,w j \x j ) 
Jrl Jr 1 : fJi 



x 





1 


) det 




/ l<j,k<L 


Uj - W k 



L" N ( Uj /2a 2 ) L» N (w k /2o- 2 ) 
L^_ 1 (^/2 ( 7 2 ) U N _ x {w k l2o 2 ) 



(3.17) 
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Proof of Eq. (2. 35) in Theorem 2.3 Since the equality (2.32) in Corollary 2.2, holds for any 
t > 0, o 2 > 0, Lemmas 3.4 and 3.6, imply that the integrand of multiple Gaussian integral 
in the RHS of (3.12) is equal to that in the RHS of (3.17). By replacing L and iV — L by n 
and N, respectively, and (ui, . . . ,ul, Wi, . . . , wl) by («i, . . . , a n , a n+ ±, . . . , a 2n ) — ol e C 2n , 
we obtain (2.35). | 



3.3 Proof of Theorem 2.4 

First we prove that (2.41) is equal to (2.33). Since 



-^w+n(aj/v / 2o 3 ) H N+n (a n+k /V2a^) 



= Hn+u-i 



«, 



H 



N+n-l 



Ctn+k 



1 



X 



1 H N+n (a j /VM)/H N+n .. 1 (a j /VM) 
1 H N+n (a n+k /V2a^)/H N+n ^ 1 (a n+k /V2a^) 



we can apply the identity (2.37) and obtain the following, 



det 

l<j,k<n 



&j O^n+k 



Hi ^ V ' 

3=1 



a, 



H N+n (aj/V2o^) H N+n (a n+k /V2o^) 

(_l)n(n-l)/2 



H 



N+n-l 



xV n ' n . . . , a n ), (a n+1 , . . . , a 2 «); 

(_]Wn-l)/2 



a n +k _ 

V2o^J S nil nLi(«n+fe - a„) 



7V+nl 



iV+n-l(7=f) 



J'=l 



if 



3=1 j 



SiN+nl ) 



where is the 2nxn rectangular matrix whose (j, fc)-element is given by a* 1 H i (aj/\ / 2a 2 ), 
1 < j < 2n, 1 < k < n. Here we use the recurrence relation of Hermite polynomials 



H n+ i(x) = 2xH n (x) - 2nH n _ 1 (x). 



(3.18) 



Then we find that 



a k ~ l H 



N+n-l 



+linear combination of jfi \{a j V / 2<r 2 ) : N + n — k<£<N + n + k — l|, 



a 



V27< 



a 



2\ (*-l)/2 

2"J 



if 



+linear combination of jif^a/v 7 ^) :iV + n- £;<£<iV + n + £;-lj, 
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for k — 2,3, ... . Therefore, by the multilinearlity of determinant, 

(TV + 



det[A N+n ^ A N+n } = a n ^ll 



t=2 



(N + n 



det[B C], 



where B and C are 2nxn rectangular matrices whose (j, /c)-elements are given by HN+n-k{&j/ V2(x 2 ) 
and H N+n+k -i(aj/V2a2), respectively. Since det [B C] = (-l) [n/2] deti< jifc < 2 n[-^Af+fe-i(tti/v / 2o 3 )], 

(_l)K2]+n(n-l)/2 = y for &ny n and 



h n (ai, a n )h n (a n+1 , a 2n ) \\ \\{a n+k - aj) = h 2n (ai, a 2n ), 

j=l k=l 

proof of the equivalence of (2.33) and (2.41) is completed. 



(3.19) 



Next we prove that (2.35) is equal to 
T(N + u + n) 1 



det 



(N + n)\ h 2n (a) l<j,k<2n 

for any v > —1. By applying the identity (2.37), we have 

1 



(3.20) 



det 

l<j,k<n 



OLj (X n + k 
_l)«("-l)/2 



n?=i Ylk=i( a n+k - a n ) 



L» N+n ( aj /2<j 2 ) L» N+n (a n+k /2a 2 ) 
Vs^iajl** 2 ) ^ N+n ^(a n+k /2a 2 ) 

det[Ajv+ n -i A N+n ], 



where is the 2n x n rectangular matrix whose (j, A;)-element is given by a k f l L»( aj /2a 2 ), 
1 < j < 2n, 1 < k < n. Here we use the recurrence relation of Laguerre polynomials 



(n + l)L" n+1 (x) = (-x + 2n + u + l)L" n (x) - (n + ^L^x). 
Then we find that 

k-l 



(3.21) 



^ k - 1 T u ( 

LX LA7 I „ li- 



ce 



N+n-l V2a2 



) = (-2a 2 ) fc - 1 ]J(N + n + v - £)L» N+n _ k (^L) 



+linear combination of ^L u t (a/2a 2 ) : N + n — k < £ < N + n + k — 1 j, 

fc— i 

(^) = (" 2(t2 )^ 1 11^ + U + t) L N + n +k -l (^) 

+linear combination of ^L u t (a/2a 2 ) : N + n — k < £ < N + n + k — 1 j, 



for k — 2,3, ... . Therefore, by the multilinearlity of determinant, 



(N + n)\ 
22 



detlW: A N+n ] = (2^)»'»-'> ( 1^1^)""' 'fj i£±!i±^det|BC], 



* r(jv + v + f) 



where B and C are 2nxn rectangular matrices whose (j, /c)-elements are given by L v N+n _ k (a.j /2a 2 ) 
and L l/ N+n+k _ l (a j /2a 2 ), respectively. Since det[B C] = (-^) [n/2] det 1 < j:k < 2 n[L u N+k _ 1 (a j /2a 2 )], 
and (3.19), the equivalence of (2.35) and (3.20) is proved for v > -1. For v e N , (2.42) 
is immediately obtained. The first expressions in (2.43) and (2.44) are obtained by setting 
v = 1/2 and v = —1/2 in (3.20) and by applying the relations (1.15) and (1.16). For the 
second expressions in (2.43) and (2.44), we use the following relations between the Hermite 
polynomials and Laguerre polynomials with v = ±1/2, 

Then the proof of Theorem 2.4 is completed. | 
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